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I. INTRODUCTION 



Strong coupling studies of lattice 2-d principal cliiral models, with the standard nearest- 
neighbour interaction 

Sl = -2NpY.^eTr [u{x) U\x+fx)] , = Wt ^ (1) 

have shown evidence of a large- iV phase transition at a finite (3c, separating the strong 
coupling and the weak coupling regions 0,0] • An analysis of the 18*^ order N = oo strong 
coupling series of the specific heat showed a second order critical behavior 

C^\(3-(3c\-\ (2) 

with the following estimates of Pc and a: Pc = 0.3058(3) and a = 0.23(3) This 
critical phenomenon is somehow effectively decoupled from the continuum limit [jS —>■ oo), 
indeed dimensionless ratios of physical quantities are reproduced with great accuracy even 
for/3</3, i|. 

A critical behavior at = oo is also present in 1-d lattice chiral models, where at 
N = oo the free energy is piecewise analytical with a third order transition between the 
strong coupling and weak coupling domains [H. In these models the parameter A^ plays a 
role analogous to the volume in ordinary systems with a finite number of degrees of freedom 
per site, and the double scaling limit describing the simultaneous approach N ^ oo and 
/9 — /9c is shown to be equivalent to finite size scaling of 2-d spin systems close to the 
criticality §,0. 

In this paper we investigate the above large-A^ critical phenomenon by Monte Carlo 
simulations, that is by extrapolating, possibly in a controlled manner, numerical results at 
sufficiently large A^, in the same spirit of the double scaling limit technique developed in the 
studies of 1-d matrix models. We performed Monte Carlo simulations of SU{N) and U{N) 
models for several large values of A^, studying the approach to A^ = oo. Some SU{N) Monte 
Carlo results at large A^ were already presented in Ref. |Q. Since SU{N) and U{N) models 
are expected to have the same large-N limit, U{N) Monte Carlo results provide further 
information and check of the N ^ oo behavior of lattice principal chiral models. 

In the continuum limit SU{N) and U{N) 2-d lattice actions should describe the same 
theory even at finite A^, in that the additional U{1) degrees of freedom of the U{N) models 
decouple. The U{N) lattice action, when restricting ourselves to its SU{N) degrees of 
freedom, represents a different regularization of the SU{N) x SU{N) chiral field theory. 
One loop calculations in perturbation theory give the following A-parameter ratios 

^ = 732 exp (0 , (3) 



where A^ and A^^ are respectively the A-parameters of the U{N) and SU{N) lattice ac- 
tions (|I|). 
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The fundamental group invariant correlation function of SU{N) models is 

G{x) = l(Tr[f/t(a:)f/(0)]). (5) 

Introducing its lattice momentum transform G{p), we define the magnetic susceptibility 
X = G{0), and the second moment correlation length 



4sin^7r/L 



^(0,0) 
G(0,1) 



(6) 



In the U{N) case we consider two Green's functions. One describes the propagation of 
SU (N) degrees of freedom: 

G(a:) = l(Tr[f/t(x)?7(0)]), 

(detf/(a;))V^ ' ^'^ 

The other describes the propagation of the U{1) degrees of freedom associated with the 
determinant of U{x): 

Ga{x) = {(det[U\x)Um)"'' ) , (8) 

From the Green's functions G{x) and Gd{x) we can define the corresponding magnetic 
susceptibilities Xd and second moment correlation lengths ^g, ^d- 

At finite N, while SU{N) lattice models do not have any singularity at finite (3, U{N) 
lattice models should undergo a phase transition, driven by the U{1) degrees of freedom 
corresponding to the determinant of U{x), and following a pattern similar to the 2-d XY 
model 1^. The mass propagating in the determinant channel Mj should vanish at a finite 
value (3d and stay zero for larger (3. Then for (3 > (3^ this sector of the theory decouples from 
the other {SU{N)) degrees of freedom, which are those determining the continuum limit of 
principal chiral models for (3 ^ oo. We recall that the 2-d XY model critical behavior is 
characterized by a sharp approach to the critical point jSxY ( the correlation length grows 
exponentially), a line of fixed point for j3 > (3xyi and a finite specific heat having a peak for 
3^13 < (3xY (see e.g. Ref. H). 



II. NUMERICAL RESULTS. 
A. The Monte Carlo algorithm. 

In our simulations we used local algorithms containing overrelaxation precedures. In the 



SU{N) case, we employed the Cabibbo-Marinari algorithm |10] to upgrade SU{N) matrices 
by updating their SU{2) subgroups, chosen randomly among the subgroups acting 

on each 2x2 submatrix. At each site the SU{2) subgroup identified by the indices i,j 
{1 < i < i < N) was updated with a probability P = j^P, so that the average number 
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of SU{2) updatings per SU{N) site variable was n = pN. In our simulations we always 
chose p ^ 1, decreasing p when increasing A^. We used p~latA^ = 9,p~2/3at 
= 15 and p ~ 1/2 at iV = 21,30. The extension to the U{N) case is easily achieved 
by updating, beside SU{2) subgroups, U{1) subgroups. In our simulations we upgraded the 
f/(l) subgroups identified by the diagonal elements of the U{N) matrix. The SU{2) and 
f/(l) updatings were performed by a mixture of over- heat-bath algorithm |Tl| (90%) and 
standard heat-bath (10%). At fixed parameter p, the number of operations per site increases 
as A^^ at large A^. 

The above algorithm experiences a critical slowing down in A^, that is keeping the corre- 
lation length fixed the autocorrelation time grows with increasing A^. This effect is partially 
compensated by a reduction of the fiuctuations of group invariant quantities when A^ grows. 
In the U{N) simulations the quantities related to the determinant channel are subjected to 
large fiuctuations, causing large errors in the measurements. 

In Tables | and y we present Monte Carlo data respectively for the U{N) and SU{N) 
simulations. Finite size systematic errors in evaluating infinite volume quantities should be 
smaller than the statistical errors of all numerical results presented in this paper. 



B. Numerical evidence of a large-N phase transition. 

Lattice chiral models have a peak in the specific heat 

^ = TvdT 

which becomes sharper and sharper with increasing A^. In Figs. ^ and |^ we plot the specific 
heat respectively for the U{N) and SU{N) models. Such a behavior of the specific heat 
should be an indication of a phase transition for A^ = oo at a finite (3c- The positions 
of the peaks Ppeak in SU{N) and U{N) converge from opposite directions, restricting the 
possible values of Pc to 0.304 ^ /?c ~ 0.309. Notice that Monte Carlo data for [3 ^ (3c — 
0.306 approach, for growing A^, the resummed 18*^^ order large-A^ strong coupling series 
of the specific heat 0; in this region, as expected by strong coupling considerations, the 
convergence of U{N) models is faster. 

A more accurate estimate of the critical coupling Pc can be obtained by using a finite A^ 
scaling Ansatz 

(3peak{N) ~ /3e + cA^"^ , (lO) 

in order to extrapolate Ppeak{N) to A^ — > cxo. The above Ansatz is suggested by the idea that 
the parameter A^ may play a role quite analogous to the volume in the ordinary systems close 



to the criticality. This idea was already exploited in the study of 1-d matrix models |T2| , 
where double scaling limit turned out to be very similar to finite size scaling in a two- 
dimensional critical phenomenon. Substituting L N and 1/u e, Eq. ([TD|) becomes 
the well-known finite size scaling relationship derived in the context of the renormalization 
group theory. Furthermore the exponent e should be the same in the U{N) and SU{N) 
models, in that it should be a critical exponent associated to the N = oo phase transition. 
Notice that the function Ppeak{N) in Eq. (^) is considered at infinite space volume. 
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In the study of ordinary critical phenomena the reweighting technique |T3|, turns out 



to be very efficient to determine quantities hke the position of the specific heat peak. In 
our work we could use this technique only for = 9, since for larger N the reweighting 
range around the point where the simulation is performed turned out to be much smaller 
than the typical j3 interval of our simulations. For > 15 (3peak{N) data and their errors 
were estimated from the specific heat data reported in the Tables | and |I|, supported by the 
direct measurements of the specific heat derivatives at each (3. 

Our estimates of pp^ak at = 9,15,21 for U{N) and N = 9,15,21,30 for SU{N) fit 
very well formula (|TUp. By a fit with four free parameters, Pc, '^u(n) ^su{n)^ found 



/9e = 0.3057(3) , 

e = 1.45(8) . (11) 

In Fig. ^ the fit result is compared with the Ppeak{N) data. A fit with two independent 
e exponents, e^^^^^ and 6^^^^^^ gave compatible results, but larger errors. Notice that this 
Monte Carlo estimate of f3c is in agreement with the determination (^ coming from strong 
coupling computations. 

We checked the finite A^ scaling Ansatz (|TD|) in the similar context of the large- A^ phase 
transition of 1-d lattice U{N) chiral models with free boundary conditions, where the critical 
point /3c and the critical exponents u and a are known: jSc = 1/2, u = 3/2 and a = —1. We 
computed the position of the specific heat peak at finite A^ finding the asymptotic behavior 
(p!0|) with e = 2/3. Details on these calculations are given in the Appendix. As already 
mentioned, from standard finite size scaling arguments the critical exponent e should be 
related to the critical exponent u: e = l/u. Notice that the critical exponents u and a 
satisfy a two-dimensional hyperscaling relation: 2z/ = 2 — a. In 1-d lattice chiral models the 
number de = 2 of effective dimensions of the large- A^ critical phenomenon is related to the 
fact that the double limit N ^ oo and /? — /3c is equivalent to the continuous limit of a 
two-dimensional gravity model with central charge c = —2. 

Since the large- A^ phase transition of the 2-d lattice chiral models is of the second order 
type, its behavior cannot be found in the classification of double scaling limits of Refs. [p!^,p|. 



which are parametrized by a central charge c < 1 implying a < 0. Moreover, unlike 1-d 
lattice chiral models, the interpretation of the large- A^ phase transition of 2-d lattice chiral 
models as an effective de = 2 ordinary critical phenomenon does not seem to be valid: 
in fact, if e = l/iy, by substituting our estimates of a and e in the hyperscaling relation 
de = {2 — a)e we would obtain de = 2.6(2). A more general thermodynamic inequality 
would give de> {2 — a)e [M. 



Monte Carlo data of x and for P Pc compare very well with the large-A^ strong 
coupling series of x (up to 15*^^ order) and (up to 14*^^ order) ||^. Fig. ^, where ^g is plotted 
versus /5, shows that data approach, with growing A^, the curve obtained by resumming the 
strong coupling series of Hi, and in particular the U{N) data, whose convergence is faster, 
are in quantitative agreement. 

Large-A^ numerical results seem to indicate that all physical quantities, such as x ^i^ud 
^G, are well behaved functions of the internal energy E even at A^ = cxd [Q. Therefore as a 
consequence of the specific heat divergence at Pc, the N = oo /3-function Pl(T) = —adT/da 
should have a non-analytical zero at Pc, that is Pl{P) ~ \P — Pc^ iu the neighbourhood of Pc- 
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By defining a new temperature proportional to the energy |T^, this singularity disappears, 



and one can find good agreement between the measured mass scale and the asymptotic 
scaling predictions in the "energy" scheme even for (3 < (3c-, where strong coupling expansion 
is expected to converge [|]. In fact strong coupling computations show asymptotic scaling 
with a surprising accuracy of few per cent |0 . 

In the U{N) Kosterlitz-Thouless phase transition driven by the determinant is 

expected at /3d > (3peak for each finite A^. Our data seem to support this picture, indeed 
after the peak of C, the magnetic susceptibility xa and the second moment correlation 
length defined from the determinant correlation function begin to grow very fast. 
In Fig. 1^ we plot Xd versus (3. Green and Samuel argued (using strong coupling and weak 
coupling arguments) that the large- iV phase transition is nothing but the large- iV limit of the 
determinant phase transition present in the U{N) lattice models PJl^. According to this 
conjecture, in the large- limit and Ppeak should both converge to Pc, and the order of the 
determinant phase transition would change from an infinite order of the Kosterlitz-Thouless 
mechanism to a second order with divergent specific heat. The large- A^ phase transition of 
the SU{N) models could then be explained by the fact that the large- A^ limit of the SU (N) 
theory should be the same as the large- A^ limit of the U{N) theory. Our numerical results 
give only a partial confirm of this scenario, we can just hint from the behavior of Xd and 
with growing A^ that the expected phase transition is moving toward /3c- The large- A^ strong 
coupling series of the mass propagating in the determinant channel has been calculated 
up to 6*^ order, indicating a critical point, determined by the zero of the Md series, slightly 
larger than our determination of f3c'- Pd{N = oo) ~ 0.324 This discrepancy could be 
explained either by the shortness of the strong coupling series of Md, or by the fact that 
such a determination of (3c relies on the absence of non-analiticity points before the strong 
coupling series of Md vanishes, and therefore a non-analiticity at Pc — 0.306 would invalidate 
all strong coupling predictions for /? > Pc- 



C. Phase distribution of the link operator. 

In 1-d principal chiral models the large-N third order phase transition is consequence of 
a compactification of the eigenvalues of the link operator 

L = U{x)U^{x + jj,) , (12) 

which are of the form A = e*^. In the weak coupling region (/? > Pc) the phase distribution 
of the eigenvalues of the link operator L, p{P,6) with 6 G (— 7r,7r], is nonvanishing only in 
the region |^| < 6c{P) < vr. The third order critical point Pc is determined by the limit 
condition 6c{P) = vr, separating the weak coupling from the strong coupling region where 
p(/5,vr)>0 §. 

In order to see if a similar phenomenon characterizes the large- A^ phase transition also in 
2-d, we have extracted from our simulations the phase distribution p(/3, 6) of the eigenvalues 
of L. Notice that p{P, 6) = p{P, —9) by symmetry, therefore in the following we will show 
p{P, 9) only in the range < ^ < vr. Large- A^ numerical results seems to support the 
compactification of the phase distribution at /3c, indeed we found p(/3, vr) ~ for /? > Ppeak 
{p{P,n) can be strictly zero only for A^ = oo). This fact is illustrated in Fig. ^, where we 
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compare the distributions p{P,0) at P = 0.300 and j3 = 0.305 for = 21, whose Ppeak — 
0.3025: the distribution values at ^ = vr (p(0.300, vr) ~ 0.010 and p(0.305, vr) ~ 0.0007) 
decrease by about a factor 15, becoming very small. Similar behaviors are observed at the 
other values of N. 

In the SU (N) models p(/5, 6) presents maxima, as Fig. ^ shows. This structure is 
absent in the U{N) models and should disappear in the large- limit, in that the height of 
the peaks with respect to the background curve should vanish. For example, the U{N) and 
SU{N) phase distributions at /? = are respectively 



and 



1 + (-1)^+^1 cos (AT^) 



(13) 



(14) 



In our SU{N) simulations we found the peak heights to decrease approximately as 1/A^. 
It is also interesting to see how the distributions p{n, /3, 6) of the generalized link operators 



Lin) 



t/(x)f/^(x + np) 



(15) 



{p{l,P,9) = p{P,9)) evolve as function of the distance n. In Fig. |^ we plot p{n,P,9) for 
N = 15, at P = 0.305 (^g — 3.79) and for various values of n. When d = n/^ ^ oo, 



p{n,f3,6) appears to tend to the (3 = distribution (14). 



D. Critical slowing down around the large- singularity. 

The large- A^ critical behavior causes a phenomenon of critical slowing down in the Monte 
Carlo simulations. At sufficiently large A^ (A^ ^ 15) and for both U{N) and SU{N) models, 
the autocorrelation times of the internal energy te and the magnetical susceptibility 
(estimated by a blocking procedure) showed a maximum around the peak of the specific heat, 
and a sharper and sharper behavior with growing A^. The increase of the autocorrelation 
times, with growing A^, was much larger around the specific heat peak than elsewhere. In 
the SU{N) simulations, te (r^) went from ~ 600 (400) at /? = 0.3025 and A^ = 21 to ~ 3000 
(2500) at (3 = 0.304 and A^ = 30 (the uncertainty on this numbers is large, they are just 
indicative). After the peak of C te and decreased, for example at A^ = 30 and /3 = 0.305 
Te — 700 and ~ 300. Similar behavior was observed in the U{N) simulations. The 
above critical slowing down phenomenon represents the most serious difficulty in getting 
numerical results around (3c at larger A^ by the Monte Carlo algorithm used in this work. At 
large correlation length increases again due the critical slowing down associated to the 
continuum limit, while te tends to be stable. 

We want to mention an attempt for a better algorithm in the U{N) case, by constructing 
a microcanonical updating involving globally the U{N) matrix instead of using its subgroups. 
A microcanonical updating of U according to the action 

A{U) = ReTi [UF] (16) 
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can be achieved by performing the reflection with respect to the U{N) matrix Umax which 
maximizes A{U): 

Unew Umax f^oici Umax i 

Notice that the determination of Umax requires the diagonahzation of the complex matrix 
F. The update does not change the action and it must be combined with ergodic 
algorithms (e.g. heat bath). We found that, at large and in the region of (3 values we 
considered, the algorithm based on the SU{2) and f/(l) subgroups performs better than 
those based on the updating (|T7|). The latter may become convenient at relatively small 
and/or for larger correlation lengths. On the other hand, at large space correlation 
lengths multigrid algorithms should eventually become more efficient, in that they should 
have smaller dynamical exponents (see Refs. [|T^,^ for some implementations of multigrid 
algorithms in the context of lattice chiral models). 



APPENDIX A 

The free energy of 1-d U {N) lattice chiral models can be written in terms of a determinant 
of modified Bessel functions 

F(iV,/5) = i^lnZ(iV,/?) = ^lndet/,_.(2iV/?). (Al) 
The specific heat can be obtained by 

The large- A^ limit of the specific heat shows the existence of a third order phase transition 
at Pc = 1/2, indeed we have 

C{oo,(3) = (3^ for (3<Pc, 

C{oo, P) = \ for p> (3^. (A3) 

The singularity at (3c can be characterized by a critical exponent a = —1. From the double 
scaling limit of the corresponding unitary matrix model the correlation length exponent 
turns out to be z/ = 3/2 |21[|. a and v satisfy a hyperscaling relationship associated to a 



two-dimensional critical phenomenon: 2ij = 2 — a. 

1-d U{N) lattice chiral models present a peak in the specific heat, whose position l3peak{N) 
should approach (3c with increasing A^. The finite A^ scaling arguments already mentioned 
in this paper lead to the Ansatz (|TU|) for the positions of the specific heat peaks. In Table ^ 
we report the values of Ppeak{N) and C{N,(3peak) as function of A^ up to A^ = 11. As shown 
in Fig. |, the large N behavior of (3peak{N) is well fitted by 

(3peak{N) = (3c + aN~' + bN-^' (A4) 

with e = 2/3, and therefore z/ = 1/e = 3/2 (a ~ 0.595 and b ~ 0.13). The result u = 3/2 
was also found in the finite A^ scaling of the partition function zeroes [|^. 
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FIGURES 



FIG. 1. Specific heat vs. /? for SU{N) models. The sohd hne represents the strong-couphng 
determination, whose estimate of the critical /3 is indicated by the vertical dashed lines. The thick 
solid lines above the peaks represent our estimates of Ppeak- 

FIG. 2. Specific heat vs. [3 for U{N) models. The solid line represents the strong-coupling 
determination, whose estimate of the critical /3 is indicated by the vertical dashed lines. The large 
solid lines above the peaks represent our estimates of /3peafe- 

FIG. 3. Ppeaki^) VS. 1/N. The dashed lines show the fit result. 

FIG. 4. vs. p. The solid line represents the strong-coupling determination, whose estimate 
of the critical /? is indicated by the vertical dashed lines. 

FIG. 5. Xd vs. /?. The vertical dashed line indicates the estimate of Pc- The filled symbols 
indicate the positions of the peak of C at A?" = 9, 15, 21. 

FIG. 6. p{P,e) for the SU{21) model at ^ = 0.300 and /3 = 0.305. 

FIG. 7. p{n, f3, 0) for the SU{15) model at (3 = 0.305 for various value of n. 

FIG. 8. Ppeaki^) vs. l/N in 1-d U{N) lattice chiral models. 
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TABLES 



TABLE L Numerical results for U{N). 



N 




L 


Stat 


E 


C 




?G 


Xd 


id 


9 


0.30 


24 


100k 


0.60374(8) 


0.284(6) 


10.50(3) 


2.058(14) 


1.64(2) 


0.6(2) 




0.31 


30 


150k 


0.56706(9) 


0.427(11) 


15.43(4) 


2.649(17) 


2.82(4) 


1.0(2) 




0.313 


30 


300k 


0.55215(10) 


0.541(11) 


18.53(5) 


2.972(14) 


3.93(5) 


1.2(2) 




0.315 


36 


150k 


0.54077(14) 


0.61(2) 


21.80(10) 


3.36(3) 


5.59(9) 


1.9(2) 




0.318 


36 


300k 


0.52128(10) 


0.66(2) 


29.47(9) 


4.08(3) 


11.3(2) 


3.0(2) 




0.3185 


36 


400k 


0.51799(11) 


0.69(2) 


31.17(11) 


4.22(2) 


12.8(2) 


3.2(2) 




0.319 


42 


300k 


0.51482(11) 


0.69(2) 


32.86(13) 


4.38(3) 


14.7(4) 


3.3(3) 




0.320 


42 


400k 


0.50816(10) 


0.66(2) 


37.2(2) 


4.73(3) 


20.5(5) 


4.3(3) 




0.323 


48 


330k 


0.49172(9) 


0.50(2) 


51.6(3) 


5.83(4) 


55(3) 


8.5(5) 




0.323 


60 


200k 


0.49166(11) 


0.52(2) 


51.7(3) 


5.79(7) 


57(3) 


7.7(7) 


15 


0.28 


18 


150k 


0.65373(4) 


0.163(3) 


6.924(7) 


1.519(4) 


1.030(8) 






0.30 


24 


200k 


0.60276(6) 


0.300(10) 


10.81(2) 


2.063(9) 


1.29(2) 






0.305 


24 


180k 


0.58405(9) 


0.396(12) 


13.09(3) 


2.346(9) 


1.57(2) 






0.308 


30 


250k 


0.56875(9) 


0.57(2) 


15.55(3) 


2.632(10) 


2.14(2) 


0.7(2) 




0.310 


30 


300k 


0.55423(12) 


0.78(3) 


18.83(5) 


2.996(11) 


3.31(5) 


1.2(2) 




0.311 


30 


500k 


0.54453(10) 


0.97(4) 


21.58(4) 


3.276(8) 


4.87(6) 


1.6(1) 




0.311 


36 


300k 


0.54470(13) 


0.97(4) 


21.52(6) 


3.276(14) 


4.85(9) 


1.5(2) 




0.312 


36 


600k 


0.53374(11) 


1.05(4) 


25.57(10) 


3.67(2) 


8.42(15) 


2.8(2) 




0.313 


36 


500k 


0.52365(12) 


0.94(4) 


30.42(10) 


4.131(15) 


14.5(5) 


3.7(3) 




0.315 


42 


200k 


0.-50920(10) 


0.50(3) 


40.0(2) 


4.95(4) 


42(5) 


7.2(7) 




0.315 


48 


300k 


0.50915(7) 


0.49(2) 


39.7(2) 


4.89(4) 


46(3) 


8.1(7) 


21 


0.28 


18 


100k 


0.65373(4) 


0.162(4) 


6.972(8) 


1.526(5) 


0.991(8) 






0.30 


24 


200k 


0.60273(6) 


0.303(9) 


10.881(13) 


2.069(6) 


1.140(8) 






0.3025 


24 


300k 


0.59390(6) 


0.361(10) 


11.869(14) 


2.185(6) 


1.220(9) 






0.305 


30 


300k 


0.58318(6) 


0.446(14) 


13.31(2) 


2.364(9) 


1.394(11) 






0.308 


30 


200k 


0.56337(15) 


0.88(6) 


16.79(5) 


2.748(13) 


2.15(3) 






0.309 


30 


300k 


0.5512(3) 


1.31(13) 


19.92(10) 


3.09(2) 


3.60(9) 


1.2(2) 




0.3095 


30 


450k 


0.5415(3) 


1.98(15) 


23.11(13) 


3.43(2) 


6.6(3) 


2.4(2) 




0.31 


36 


300k 


0.5337(2) 


1.28(10) 


26.23(11) 


3.75(2) 


11.1(5) 


3.4(3) 
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TABLE II. Numerical results for SU{N). When more than one lattice size appears, the corre- 
sponding results were obtained collecting data of simulations at the reported lattice sizes (which 



were, in all cases, 


in agreement within the 


errors) . 
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30 
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24 
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36 
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u.ouz 


36 
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U.OOOUU^O J 
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O. ±0±yiJ ) 


0.305 


36 
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45 
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24 
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0.3025 


30 
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17.02(4) 


2.784(8) 


0.3025 


36 


500k 


0.56491(11) 


1.04(5) 


17.11(4) 


2.800(10) 


0.303 


36 


500k 


0.55959(9) 


0.96(4) 


18.38(3) 


2.936(9) 


0.305 


30 


500k 


0.54100(8) 


0.72(2) 


24.14(5) 


3.526(8) 


0.310 


42 


240k 


0.51548(6) 


0.41(2) 


36.66(12) 


4.61(2) 


30 0.300 


24 


150k 


0.59927(8) 


0.38(2) 


11.35(2) 


2.114(7) 


0.3025 


30 


200k 


0.58479(10) 


0.79(5) 


13.24(3) 


2.338(7) 


0.303 


30 


200k 


0.58007(15) 


1.00(8) 


13.99(4) 


2.433(10) 


0.304 


24 


500k 


0.5625(4) 


2.4(3) 


17.55(7) 


2.857(10) 


0.304 


30 


500k 


0.5632(3) 


2.3(2) 


17.40(7) 


2.829(8) 


0.305 


30 


200k 


0.5466(2) 


1.05(10) 


22.13(5) 


3.320(12) 





TABLE III. Ppeak{N) and C{N,l3peak) versus N for 1-d C/(A^) lattice 


chiral models. 


N 


f^peak 


C{l3peak) 


2 


0.930889 


0.29461215 


3 


0.818356 


0.27992604 


4 


0.758001 


0.27269388 


5 


0.719664 


0.26839003 


6 


0.692846 


0.26553250 


7 


0.672876 


0.26349442 


8 


0.657337 


0.26196545 


9 


0.644848 


0.26077452 


10 


0.634554 


0.25981956 


11 


0.625899 


0.25903594 
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Figure 6 



P(p,e) 



0.30 



0.25 



0.20 



0.15 



0.10 



0.05 



0.00 



0.0 



1.0 



2.0 



P=0-300 
P=0-305 




3.0 



e 



Figure 7 




Figure 8 




